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Abstract 

We discuss hypersurface motions in Riemannian manifolds whose normal velocity is 
a function of the induced hypersurface volume element and derive a second order 
partial differential equation for the corresponding time function t{x) at which the 
hypersurface passes the point x. Equivalently, these motions may be described in 
a Hamiltonian formulation as the singlet sector of certain diffeomorphism invariant 
field theories. At least in some (infinite class of) cases, which could be viewed as a 
large- volume limit of Euclidean M-branes moving in an arbitrary M + 1-dimensional 
Riemannian manifold, the models are integrable: In the time-function formulation the 
equation becomes linear (with t{x) a harmonic function on the embedding Riemannian 
manifold). We explicitly compute solutions to the large volume limit of Euclidean 
membrane dynamics in IR'' by methods used in electrostatics and point an additional 
gradient flow structure in IR". In the Hamiltonian formulation we discover infinitely 
many hierarchies of integrable, multidimensional, A''-component theories possessing 
infinitely many diffeomorphism invariant, Poisson commuting, conserved charges. 
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1. In physics, and mathematics, surface motions are usually considered in- 
dependent of the parametrisation, the most prominent example perhaps being 
the 'flow by mean curvature' which has been of intense mathematical interest [Q] 
as well as of physical importance . When studying the dynamics of relativis- 
tic membranes, on the other hand, one encounters [Q in a partially gauge-fixed 
formulation, a very interesting equation (describing the time-evolution of a 2- 
dimensional surface in IR^) which is parametrisation- depenrfent (meaning that 
different parametrisations of the initial surface Eq generically lead to geometri- 
cally different shapes at later times, t (in the case at hand there is a residual 
symmetry group of area-preserving diffeomorphisms, but any reparametrisation 
whose Jacobi determinant is different from 1 will generically change the motion); 
alternatively, the equation (though first order in time) may be viewed as requir- 
ing the initial unparamctrizcd shape Sq and the initial unparametrized normal 
velocity field on Eg as initial conditions (the parametrisation of Eq, and of Et, 
then follows, up to area-preserving diffeomorphisms of Eq -that don't change 
the geometry of the motion- from the equation). Moreover, it was shown in 
that this seemingly parametrisation dependent equation can be obtained from 
the singlet sector of a diffeomorphism invariant Hamiltonian field theory. In 1^ , 
on the other hand, some surface motions in IR'^ were shown to be best described 
by a (second order partial differential) equation for t(x), the time at which Ej 
contains x. 

In this letter, we significantly generalize this situation (to higher dimen- 
sions, curved embedding spaces, and more general dynamics) by considering 
a rather general class of parametrisation dependent hypersurface motions in 
Riemannian manifolds. For these hypersurface motions we shall give an equiva- 
lent diffeomorphism-invariant (constrained) Hamiltonian formulation as well as 
derive a second order partial differential equation for the time- function t(x), 
which in some cases (that may be viewed as large-volume limits of Euclidean 
A/-branes or, equivalently, as highly generalized, multilinear, field theoretic ex- 
tensions of Nahm's matrix equations [^) turns out to simply be At = where 
A is the Laplacean of the embedding Riemannian manifold in which the motion 
takes place. In the particular case of three-dimensional Euclidean space the 
solution of the Laplace equation for the time function becomes easy to handle 
thanks to the formal equivalence with the potential equation of electrostatics: 
we compute explicit solutions of these surface motions determined by singular 
one-dimensional membranes at t = —oo and t = +oo (corresponding to neg- 
atively and positively charged composed loops of wire, respectively) between 
which two-dimensional surfaces develop in time. 

Finally, and perhaps most importantly, we are led to an infinite set of hi- 
erarchies of multidimensional integrable diffeomorphism invariant iV-component 
field theories for which we are able to provide a complete set of Poisson-comniuting, 
diffeomorphism invariant, conserved charges. 

2. Let E be an orientable compact connected manifold of dimension M with 
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a fixed volume form p. We shall denote co-ordinates on E by (^'", 1 < r < M. 
Let (A/", C) be an orientable Riemannian manifold of dimension N = AI + 1 on 
which we shall denote co-ordinates byx*,l<i<M-|-l and whose volume form 
induced by C will be denoted by loc^. Consider the set of smooth immersions 
a; : E ^ A/", inini{I] , Af) . For each immersion x there is a unique outward 
normal field n[x] : S — > TAf (depending on x and its first derivatives) attached 
to the parametrized hypersurface x : S ^ J\f. Moreover, let x*( denote the 
Riemannian metric on S obtained by pulling back ^ to E by x and denote by 
LOx'^ the corresponding volume form on E. It follows that there is a unique 
smooth positive function y/ g[x]/ p : E ^ IR depending on x, its first derivatives, 
and the reference volume form p such that Wx*c — {\/ 9[^]/ p)P- Finally, let a 
be a positive smooth fuction on (a suitable open interval of) the positive real 
line with nowhere vanishing derivative. Consider now the hypersurface motion 
(— e, e) imm{T,,Af) : t ^ x{t) in N described by the following differential 
equation: 

dx I 

"(v5N/p)nN . (1) 



dt 

In local co-ordinates this equation looks as follows: writing 



d^}"^ dip 

for the induced metric x*C, we get 



dx^ dx^ , , , 

QAx) (2) 



= a{^y det{grs) / p) 



1 r)T*i fjq-iM 

C^(x)yd^^(aO(^)e.n---.M^'-^---'-'-^7r:- • • • TT— (3) 



The following choices of a are of particular interest 



Lorentzian A/— brane : x = y I — {\/ g[x]/ p)'^ n[x] (4) 
Euclidean A/— brane : 



\/i79^/pf~^ n[x] (5) 



and the large- volume limit of (^, 

i ^ i\/9M/p) n[x]. (6) 
The small- volume limit of (^, 

X = n[x] (7) 

(which may be called the Optical model) is not contained in this class (since the 
prefactor a of the normal is constant). Nevertheless, it can be solved directly 
since (0) is easily seen to imply the following equation of free motion: 

- , , dx^ dx'^ , . 
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(where F^j, denote the ChristofFel symbols of () which can be solved as soon as 
the geodesic flow of {J\f, C) is exphcitly known (e.g. for flat IR^, the iV-sphere, 
etc.). Some of the following results will also remain true for (|^). 

As in we are trying to derive a differential equation for the time-function: 

Theorem 0.1 Let E, p, J\f, C o.nd a he defined as above. 

1. Suppose that there is a positive real number e and a solution x of the above 
equation (Q) such that the map x : (— e, e) x S — > A/" : (i, i-^- x{t, (p) is 
a diffeomorphism onto its image, AC, which is an open neighbourhood of 
the initial hypersurface x(0, S). 

Then the first component of the inverse map x~^ : Afe (— e, e) x E 
which we shall call the time function r satisfies the following equations 
where Vr := (^'^{dr) (V*t = Q^djr) denotes the gradient of t, \v\ := 
y/ <^{v, v) = yfCjj^fHfl for any tangent vector v to M, and the symbol ";" 
stands for covariant derivative. 

Vr 

M (9) 



iVrl 



|Vr|(x) = , 1^, ^ , (10) 



= iVrpCV.y. + (/3(|VT|)-l)T;,yVr'VT-''. (11) 
where (3 is defined by: 

m ■■= -z§-^ia-'il/z)). (12) 

Eqs d) and ^ remain true for positive a whose derivative may vanish. 

Conversely, suppose that the smooth function t : Af IR is a solution 
of the second order partial differential equation ( pl[ ) and that there is a 
positive real number e such that all its level surfaces t = c for c G (— e, e) 
are diffeomorphic to E. 

Then there is a parametrisation of the zero level surface Eq of t , Xq : T, 
Eo and a positive real number r such that eqn ( jJ^ is satisfied for xq and 
rp. Moreover, let^t denote the flow of the vector field X :— Vt/|VtP, i.e. 
the solution of the ordinary differential equation d^t{x) / dt — X{<^t{x)) 
with initial condition ^q{x) = x. Then for each parametrisation xq of 
the zero level surface which satisfies the above condition ( [7^ the map 
x : (— e, e) X E — > A/" defined by 

x{t,^):=^t{xQ{v)) (13) 
satisfies equation ( with p replaced by rp ) with initial condition . 
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Proof: 1. It is obvious that the gradient of r wiU be orthogonal on the surfaces of 
constant time whence the gradient of r is proportional to the surface normal 
n[x] which gives (^. Using the inverse function theorem 

_ dx^ dr dx^ dip'' 
dt dx^ difi^ dxi ' 

contracting with (^ikn[x]'' , and using equation (|l]) and the definition of the surface 
normal we get the second equation ([L0|). This leads to 

Vr 

• ^ (x) (14) 



In order to obtain a differential equation for the time function r in the case where 
the derivative of the function a does never vanish we differentiate equation ( p^ 
with respect to time and use equation to get rid of x. Thereby the time 
derivative of the left hand side of can solely be expressed by second covariant 
derivatives of the time function. To do the same for the right hand side we use 
the equation 

obtain by the chain rule 

drX' = d,{VT/\VTf)drX' 

observe that the orthogonal projection n onto the tangent spaces of the hyper- 
surfaces can be expressed in two ways by 

g^'^drx'dsX-' — TT^-' (x) = C'^(^) ~ n[a;]'n[a;]-' 

(see e.g. Q), and replace the normal vectors by the normalized gradients of the 
time- function ^ which gives the second order equation (pT|). 

2. Denote by Eo the level surface r = and by i : So ^ A/" the canonical 
inclusion. Let tJi*^ denote the volume form on Eq induced by the metric (. 
Consider the modified volume form 

1 



on Eo. Let r be the unique positive real number such that the two following 
integrals are equal: 

1 



By Moser's lemma (see [|ll|) there is a diffeomorphism xq : E Eo such that 
the two volume forms are diffeomorphic, i.e.: 

MVTl(i)-' y\VT\{xo)' 

whence eqn (^o|) is satisfied for t — and rp. 
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Observe now that for all points y £ Eo we have 

r($i(y)) =t 

which can be seen by differentiating the left hand side with respect to t and 
showing it to be equal to 1. When we insert y — xo{ip) in this equation and 
differentiate with respect to i/p'" we see that the gradient of the time function is 
always orthogonal to the surfaces $t(a;o(S)). Hence Vr(a;(t, <^))/| Vr(a;(f, ^))| 
is equal to the hypersurface normal n[xt]{ip)) whence we have the differential 
equation 

., , VT(x(t,Lp)) 1 r 1/ \ 

We have already shown above that the function (for fixed <^ £ S) 
fit) ~ ai^gr^t, ^)/{rp{^))) \V{x{t,ip))\ 

equals 1 for t = by construction of a^o. Upon differentiating / with respect to 
t, eliminating x by the equation above and using the second order equation (^l|) 
we get the first order time dependent differential equation for /: 

. _ r;.jVr'VT^/3(|Vr|) , 1 

lVr|* ^""^'""''hiWrDixit,^)) pO-Ei1ML^))^^'> 

Since the right hand side smoothly depends on / in a neighbourhood of / = 1 
the obvious constant solution f{t) — 1 is unique. This proves eqn (^) for all t. 

Q.E.D. 

In other words, the open neighbourhood AC will be foliated by the level 
sets of r which are nothing else but the surfaces of constant time. Eqs (^) and 
( p^ have also been considered in and and only use the fact that x is 
proportional to the unit normal. 

Note that for nonconstant a the above surface motion is not reparametrisa- 
tion invariant, which is here reflected in the fact that the time function satisfies 
a second order partial differential equation whose normal derivative at the zero 
level surface encodes the "volume part" of the parametrisation, i.e. the induced 
volume form compared to a reference volume. On the other hand, for the Op- 
tical model (0) the function a is equal to one: this model of surface motion is 
invariant under reparametrisation, and the differential equation for the time- 
function is first order, viz. eqn (^) for a = 1 which is the eikonal equation 
|Vrp — 1 known in optics. 

The second order equations for all the time- functions dealt with in Theorem 
1| can be derived from a Lagrangean: define the real-valued smooth function 
F as any solution of the first order differential equation: 
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where C is an arbitrary nonzero real number. Define the Lagrangean £ for the 
time-function r as: 

C{dT) := ^det(C.,) F{\Vt\) (16) 

A straight forward calculation shows that the Euler Lagrange equations produce 
(0) for aU F solving 

For the model we can choose the Lagrangean £(9t) equal to -ydet^^^j Vr|' 
and get the Laplace equation 

Ar = (17) 



where A is the Laplacean operator —. — 1 9; ( -v/det ( Qj ) dj . This means that 

V dct(Cij) 

-despite its complicated nonlinear appeal- the model (^ is integrable at least for 
all those Riemannian manifolds for which the Laplace equation can be solved 
in a sufficiently explicit way, e.g. for IR", the sphere S*" and hyperbolic space 
H". We shall discuss the surface motion in IR"^ (for which a linearizability was 
already noted in without using the time function) according to its harmonic 
time-function in more detail in the next section. 

In numerical relativity theory the foliation of a given spacetime into spacelike 
hypersurface which are level surfaces of a (Lorentz) harmonic time-function has 
been dealt with by C. Bona and J. Masso in This method of 'harmonic 
slicing' avoids the occurrence of spacetime singularities on these surfaces. 

For the Lorentzian membrane we can choose the Lagrangean equal to the 
function •y^det(Cy)^|Vrp — 1 and get the Euler-Lagrange equations 

(1 - |VTp)AT-|-r:yVT*Vr^' = (18) 

whereas for the Euclidean membrane we can choose C{dT) equal to the function 
A/det(Cij)-\/l + iVrp and get the equation 

(1 + |Vt|2)At - r;y Vr'Vr^' = 0. (19) 

Note that the global description of surface motion by parametrisation may 
drastically differ from the description by the level sets of a time function: al- 
though being equivalent for short time intervals (and some technical assump- 



tions, see Thm 3.1) the former allows the membrane to be at the same position 
at different times, but fixes the topological type, whereas the latter allows for 
varying topological type, but every point in space is contained in at most one 
level surface. In both pictures one has to admit violations of the regularity of 
the mappings (i.e. points where the Jacobians cease to have maximal rank) in 
order to produce interesting situations. 

3. In this section we should like to discuss some solutions of the model (||) 
in IR'^ in terms of its harmonic time-function. A comparison with the Euclidean 
membrane model (^ shows that (|^) is an approximation of the Euclidean mem- 
brane model in the regime where the induced spatial membrane area element 
^Jdet{grs) is large compared to a reference area element p. 
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The Laplace equation (|T^) for flat IR'^ can be treated by methods known in 
electrostatics (see e.g. Q). For example, consider a piece of wire of length L 
which is uniformally charged with line charge \ = q/ L and which is lying along 
the 2-axis between —p := L/2 and p: its electrostatic potential r is given by 

r (f ) = A log ( V^^^^^±1±P) (20) 

^ -\- 1p + {z — p)"^ + z — p 

This function is harmonic outside the interval [—p, p] along the z-axis on which 
it becomes ±00 according to the sign of A and tends to zero for large distances 
from the origin. It is known that all the equipotential surfaces of this function 
are axially symmetric ellipsoids whose focal points are the end-points of the 
piece of wire. In case A is negative we get the -at least in electrostatics- unusual 
picture of the equipotential surfaces as surfaces of constant time: for t = —00 it 
is a singular line which blows up into bigger and bigger ellipsoids, and at r = 
the level surface is the sphere at spatial infinity. 
This situation can be generalized a bit: let 



:5^^IR^ l<i<m, l<j<n (21) 

be closed curves in IR^ which are either points or immersed with only a finite 
number of self- intersection points. Suppose furthermore that any two different 
curves do not intersect. Moreover, let A~,A^ : ^ IR be positive smooth 
functions where 1 < z < m, ^ < j < n. On IR'^ minus all the images of 
the closed curves, a?]" (S*^), . . . , x^(5'^), ^^(5^), . . . , a;+(5'^) we can define the 
following harmonic function r: 



i=l 



(It is harmonic because x 1— > j^z^ is harmonic for all x ^ y.) From the point 
of view of electrostatics this a finite collection of closed curves carrying either 
positive or negative line charges. Due to the presence of the denominators in 
the above formula for r it is clear that there is a negative real number —t^ of 
very large absolute value such that the equipotential surface t = —to breaks up 
into m connected surfaces located each in the vicinity of one of the "incoming 
curves" x~(S^\ The following intuitive argument shows that each connected 
piece is diffeomorphic to a Riemann surface of genus equal to the number of loops 
which are generated by the self-intersections. In other words: the level surface 
piece near such a curve looks like a two-sphere if the curve is just a point, like 
a two-torus if the curve is immersed with no self-intersection, like a surface of 
genus two if the curve looks like the figure 8, and so forth. Indeed, since x~ {S^) 
is compact we can cover it by a finite number of balls of diameter L where L 
can be chosen arbitrarily small. Now suppose that the map x~ is well-behaved 
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enough so that we can approximate it by a sequence of straight pieces of wire 
each carrying a constant hne charge equal to the mean value of A~ over the 
piece of the curve. Formula (|2^) gives an approximate solution of the potential 
T. From this formula one can infer that if the distance of the point x to the 
curve is very small compared to L the potential is approximately proportional 
to the natural logarithm of the distance. Hence for "very early" values of t 
the connected components of the equipotential surface are approximated by 
surfaces of constant distance from a given curve which visibly have the required 
topology. One has the same argumentation for very large or "very late" values 
of r: here the connected components of the equipotential surfaces are certain 
Riemann surfaces centered around the "outgoing curves" x^{S^). 

We can interpret this geometrical picture as a theory for propagating mem- 
branes with prescribed change of topology: if the time equals — oo one has m 
incoming membranes shrunk to one-dimensional objects whose homotopy and 
initial "shape" are encoded in the curves sc~ with line charge functions A~. 
These m membranes will grow bigger and possibly melt together (which can be 
interpreted as interaction) and for large positive time values the surface of equal 
time decomposes again into n connected components which will eventually cen- 
ter around the outgoing curves with final "shape" described by and will 
be equal to these curves for t = +00. 

In order to get more general solutions one has to solve the Dirichlet problem 
for the potential r outside m conducting Riemann surfaces held at value —to 
and n conducting Riemann surfaces held at value +to for a positive real number 
to whose solution in principle exists and is unique. 

In order to get explicit solutions for the model (^ in higher dimension- 
sional Euclidean space IR*^^^ we can easily generalize the above construction 
(p2[): we have to replace the curves xf{(p) and by (at most) M — 1 dimensional 
surfaces xf{ip^, . . . , (p^'^~^) equipped with higherdimensional "line charge" func- 
tions and the line integral over l/\x — x~{ip)\ by a surface integral (of codimen- 
sion two) over l/\x — x~{(p^, . . . ,ip^^~^)\^'^~^. More generally, one has to ex- 
plicitly solve the higher dimensional analogue of the above-mentioned Dirichlet 
problem whose solution in principle exists and is unique. 

Finally one should point out that (ph may also be written as 




(23) 



with 



(24) 



W[x] 



M + 1 



1 



■III •••ijf 



e 



(25) 
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in ]R", i.e. are gradient flows with respect to the volume functional W and the 
metric H^^ {ip, (p). 

4. Consider now a Hamiltonian field theory with fields from E into the 

cotangent bundle T*M of M where the conjugate momentum field p is supposed 
to be densitized. ^ Let the Hamiltonian be of the form 



H{x,p) := / MAK^,) (26) 



where /i is a smooth function of one variable, ^/g :— ^/g[x] :— \J det g^-s is defined 
as in (|2|), p^ :— piPjQ-^ {x)^ \p\ :— \fp^, and u will be short for -|^. |^ 
From the canonical Poisson structure 

. , [ M i SG SF ^ 

^ '^"7/ '^^feH<p)<5p.M Sx^{^)Sp.{v)' ^^^^ 

we obtain Hamiltonian equations of motion 

Cjj (x)9^((/i - uh')^g''''dsX^) 
^C.^j{x){h - uh')y^g'''drx''dsx'r{i{x) 
-ld.e\x)h'^ (29) 

Cr{x,p) :— PidrX^ (30) 



dt 



Defining 



(the infinitesimal generators of diffeomorphisms) it is easy to check that ( p8[ ) 
und (E9h imply 

(TJ is invariant under diffeomorphisms : E ^ S), as well as 
B hh' 

^(V?/.) = a.((— - h'')g^^C,ix,p)). (32) 



^Although we shall be always working in canonical x' and pi co-ordinates (a bundle chart) 
on T*Af this has a global meaning: the pair {x,p) is a smooth vector bundle homomorphism 
A*^TS — > T*M over the map x : E — > A/', i.e. xoA'^^tj; = t^o {x,p) with the obvious bundle 
projections. 

''in order to avoid clumsy notation we shall henceforth suppress the arguments (/3 £ S in 
the integrals, the field x and its derivatives in ^[x] and n[x], and the argument u of h and 
h' from time to time. 
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For solutions {x{t),p{t)) of the equations of motion satisfying Cr{x{0),p{0)) — 
the r.h.s. of cqn ( |3^ ) vanishes identically (for all t due to (|3l])), hence (for such 
solutions) 

^[xit)]h{^^) =: p (33) 

with p some time independent density on S. Inverting (|3^) to obtain as a 
function of and noting that for — Af = 1 the condition 

Cr{x,p) = p.drx^ = 0, 1 < r < M (34) 

for X an immersion means that Q-'pj must lie in the direction normal to the 
hypersurface defined hy Lp x{t, ip), one sees that ( ^8|) can be rewritten in the 

m 



form mentioned in section 2 (sec (|l 



dt ' p 

with the functions a and h being related via 



- a{^)n (35) 



a(z)z25,(/i-i(i)) = -l. (36) 



man 



z 

In particular, for h{z) = \f2z we get oi{z) — z ^ and 

/ d*V (4.9Np2)i (37) 
For h[z) — Vl + -2^ which corresponds to q;(z) = \/l — we get the Hamilto- 

H[x,p] = [ d^V Vp' + (38) 

of the relativistic Af-branc in Af + 1 dimensional Minkowski space. h{z) = z 
implying 01(2;) = 1 gives the Hamiltonian 

H[x,p] = J d^^Lp (39) 

of the Optical model (Q). 

Let us for a moment concentrate on ( pZh: Motivated by the fact that the 
equation for the time function is linear (see |l7|) and the observation made 
for TV = 3, that a Lax pair formulation of dx/dt = -y/gn implies the time- 
independance of / Q{x{t, (p^, (p'^)dip^d(p'^ for harmonic polynomials Q{x^ ,x^ , x^), 
one is led to the conjecture that for all = Af + 1, p, and Q{x) any harmonic 
function of x, 

Q[x] := / P Qix) (40) 
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will be independent of time if the hypersurface St = a;t(S) evolves according to 
(|). The proof is a trivial application of Gauss' Theorem: assuming that for all 
t the moving surface St is a boundary of some open set Vt of Af and writing 
dSl for the surface element we get 

dCl[x{t)] f /.^^^^' 

- ' d'^'ipp diQ[x{t))- 



dt .u^ " ' " at 



dSl diQ{x{t)) 

St 



d^a;v/det(Cfci) d,Q) 
f d^'x^detiCki) AQ = 0. 



In the Hamiltonian formulation ( |37| ) the nondynamical density p{(p) is simply 
replaced by the Hamiltonian density 

i.e. TClXjp] — {g[x]p^)i in the case of (|37|). Indeed, for any harmonic Q the 
following functionals on phase space, 

Q[x,p] := / d^V {9p')-^Q{x), (42) 
Js 

will be time-independent for solutions {x{t),p{t)) satisfying (p^), (p9|), and the 
constraint (|34|). 

Do the quantities ( |4^ ) Poisson-commutc (on the reduced space)? The an- 
swer is 'yes' (in fact they form an infinite, presumably complete set of Poisson- 
commuting, reparametrisation invariant charges, thus providing the solution of 
a non-trivial diffeomorphism invariant field theory)"but let us first make some 
observations valid for all theories of the general form (p6[), i.e. arbitrary h, 
N, M, and ^: With the canonical Poisson structure ( p7| ) one has the following 
Poisson brackets where /, / are arbitrary smooth real- valued functions on S and 
X, Y are arbitrary vector fields on S: 

{ / ^""Cr, / d^^if yC,} = / [X,YYCr (43) 
je Jh js 
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(44) 

{ / m, I d^v y'Cs} = - I y^dsfu (45) 

where denotes the Lie bracket dsV^X" ~ dsX^'Y'' of X and F. 

While (^), together with (plj), imply that in principle the Hamiltonian 
theory ( |26|) may be reduced to a Hamiltonian theory on a difFeomorphism- 
invariant phase space (involving only N — M fields and N — M momenta) one 
should note that on the reduced phase space the density Ti ( ^ ) can no longer be 
defined, unless integrated over particular functions. This explains that although 
?i = ( ^2|) for solutions of ( |28| ) and (^^ satisfying the constraint ( |3l| ) , the Dirac 
bracket of H, with H (see e.g. for the case (^8|)) would not give zero, due to 
(^5|). In particular, 7i is not a conserved quantity in the Hamiltonian theory. 

The charges Q (|4|), however, are conserved on the reduced phase space: 
Q[x,p] — J^d^'^ (p 'HQ{x) commutes with Cr, 

{ I d^V 'HQ[x), j yCr) = - I Y''drQ{x)n 

d*V nd,Q{x)Y''drx' = (46) 



E 



(irrespective of Q, /i, and iV — M; just use ( ]4^ ) with / = Q{x) for the first term 
on the r.h.s of (p6[)), as well as with H (provided Cr ^ 0, N = M + 1, and 
hh' = const). Moreover, 



{ / nQix), f d^'ip HQ{x)} 
je Je 



hh' _ _ 

( h'^)g'-'{Q{x)dsQ{x) - Qix)dsQ{x))Cr 

E U 

d^VVsN /»'/if^C^(:r)(Q(:r)9,Q(:r) - g(a;)a, Q(a;)) 

E \P\ 
M 



const. / d ^p^/g[x] n^[x][Q{x)djQ{x) — Q{x)djQ{x)) 

JE 

const. / d^xy^i d, [QdjQ - QdjQ)) = (47) 



where again we have assumed that Sf is the boundary of Vt. In (47), the first 
equality is general (using (^), with / Q{x), f = Q{x), for the first term), the 
second equality holds if = (1 < r < A/) and TV = M + 1 (so that C^pj is 
proportional to the surface normal), and the last equality holds if Q and Q are 
harmonic functions (it actually holds as long as QAQ — QAQ — 0, which is a 
much weaker requirement). 

The preceding observation has far-reaching consequences. It not only proves 
the Poisson-commutativity of the conserved charges Q given in (B3) (-hence the 
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integrability of (|37|); the freedom of choosing Q{x) inJR^ e.g. as J2hLo \m\<i limYimi 
with Yim{x) being the sohd spherical harmonics, neatly matches the degrees of 
freedom to be expected from the one single field that is left over after the Hamil- 
tonian reduction-) but -confirming that integrable theories come in hierarchies- 
shows that any of the charges (^2|) with AQ — may be considered as the 
Hamiltonian of another diffeomorphism invariant field theory (possessing ( p2[ ) 
as commuting conserved charges). 
Due to 

'"d*V Q{x){'ig[x]p^)i = I d'^'ip iig[x]p^)i =: H[x,p] (48) 

corresponding to Q{x) — 1 for a Riemannian manifold (A/", C) with conformally 
equivalent metric 

4(x) :=(g(x))^C.,(a:) (49) 

(implying Vff = iQ{x)f^^, f = p.pjC^ = PiPjC^Q'^] N ^ 2] the very 
special case N — 2 can be dealt with separately), and the fact that under the 
conformal change ( ^9|) the scalar curvature changes according to 

= Q-^R (51) 

one sees that each hierarchy (of integrable Hamiltonian systems, resp. harmonic 
time-functions) consists of hypersurface motions in Riemannian manifolds hav- 
ing conformally equivalent metrics and satisfying 

RC~,=RQj. (52) 

This intriguing observation should have many important consequences. Taking 
e.g. Af = IR^ and = Sij the above suggests that the solutions of the linear, 
but nontrivial, equation 

V • (g2(f)VT(f)) = (53) 
resp. the corresponding nonlinear equations 



are related to 



resp. 



-=Q(.) n (54) 



VV(f) = (55) 
dx \/det(aH?~aIf) ^ . 

m ^ -p " 
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via some generalized Backlund transformation for all Q{x) satisfying 

V2g(f) = 0. (57) 
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